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Abstract 

Five different versions of the three-dimensional (3D) reduction of the Bethe- 
Salpeter (BS) equation in the instantaneous approximation for kernel of BS 
equation for the two-fermion systems are formulated. The normalization con- 
dition for the bound-state wave function in all versions are derived. Further, 
the 3D reduction of BS equation without instantaneous approximation for the 
kernel of BS equation is formulated in the quasi-potential approach. Except 
of the Salpeter version, other four versions have the correct one-body limit 
(Dirac equation) when mass of one of constituent fermions tends to infinity. 
Application of these versions for investigation of the different properties of 
the qq bound systems are considered. 

PACS numbers: 03.65.Ge,03.65.Pm,11.10.St,12.39.Ki,12.39.Pn 
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I. INTRODUCTION 



After having firmly established the quark structure of mesons and baryons, there nat- 
urally arises the question: how to describe the properties of hadrons in terms of explicit 
quark and gluon degrees of freedom. The main feature of QCD at low energy - confine- 
ment of quarks and gluons into a colorless bound states - is still understood very little. For 
this reason, one has to resort to various kinds of QCD-inspired models. The simplest one 
is the so-called constituent quark model, where quarks have a given "constituent" mass, 
and the interactions between the "constituent" quarks within mesons qq and baryons qqq 
are described by "confining potentials" , growing to infinity at the infinite quark separation. 
At the first stage, this intuitive picture has been implemented within the non-relativistic 
approach. Despite the evident success of the non-relativistic potential model, it has been 
understood long time ago, that one has to include the relativistic effects, at least when 
describing hadrons consisting of light u, d, s quarks. Field-theoretical Bethe-Salpeter (BS) 
equation provides a natural basis for a relativistic generalization of the potential model, 
where both light u, d, s and the heavy quark c, b bound states can be treated on the equal 
footing. A more sophisticated approach is based on a coupled set of Dyson-Schwinger (DS) 
and BS equations, that can be derived at QCD level [0-0. In such an approach, one uses 
a model gluon propagator that through the solution of the DS equation leads to the quark 
propagator which is an entire function in a complex p 2 plane and therefore is believed to 
correspond to the confined quark. A full content of underlying QCD symmetries which are 
important at low energy, can be consistently embedded within this approach. In particular, 
the Goldstone bosons are properly described, and in the limit of the vanishing quark masses, 
the masses of Goldstone bosons obtained through the solution of the coupled DS and BS 
equations, also vanish (Note that it is not the case in the simple potential-type models with 
quarks having the constant "constituent" mass). 

In the following, we shall review the potential model based solely on the BS equation, 
which is the subject of intensive investigations during last twenty years. 

II. BETHE-SALPETER EQUATION FOR THE TWO-FERMION BOUND STATE 

To set up the notation, in this section we give a brief survey of the covariant BS approach 
to the two-fermion (fermion-antifermion) bound states. In order to derive the BS equation 
for the two-fermion bound state, we consider the full 4-fermion Green function G which in 
the momentum space is given by 

G{pi,p 2 ;p' 1 ,p' 2 ) =i 2 J dx^dx^dx'z ^1+^2-^-^4 x 

x (01^1(2:1)^(^)^1(^1)^2(4)10) > (2- 1 ) 

where, for simplicity, the fermions 1 and 2 are assumed to be distinguishable, and the spinor 
indices are suppressed. 

The Green function satisfies the BS equation in the momentum space 

G = G + G KG = G + GKG , (2.2) 
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where Go stands for the free 4-fermion Green function (the direct product of two fermion 
propagators), and K denotes the kernel of BS equation, given by the sum of all two-particle 
irreducible Feynman graphs. 

In the momentum space, it is convenient to define the center-of-mass (cm.) and relative 
4-momenta according to the following relations (with arbitrary a and 

P = Pi + P2, P = Ppi - ocp 2 , a + (3 = 1 , 
or pi — aP + p, p 2 = f3P — p. (2.3) 

For the basis vectors in the momentum space, the following notation is used. 

\Pi) ® |P2> = \P1P2) = \Pp) = \P) ® \p) ■ (2.4) 
These vectors satisfy the completeness and orthonormality conditions 

/W^W = l fcri = l,2, /|P)g(P| = l, /irf^OH-!, (2.5) 

(p,lp;>=%(2ir)V(p, - P ;), = (2vr) 4 « 4 (P - P') , (p|p'> = (2ir)W(p - p') . (2.6) 

In these notations, we can write 

(Pp\0\P'p') = (2tt)V(P - P') (p\0(P)\p') = 0(P;p,p') 

Further, 

(p\G (P)\p') = G (P;p,p') = (2nf5\p-p')G (P;p), (2.8) 

G (P;p) = Si( Pl ) ® S 2 (P2) = -{A + m 1 )®{p' 2 + m 2 ) g (P;p), (2.9) 

where Si(pi) = — m,) -1 stands for the free fermion propagator with the mass rrii, and 
the quantity g (P;p) is defined as follows 



, Q = G, G , K. (2.7) 



g (P;p) p2 _ m2 + iQp2 _ m2 + t0 p 2 o _ w 2 + t0p 2 o _ w 2 + t0 > ( 2 - 10 ) 



with Wi = Jm? + pf. 

The stable bound state with the mass M B in quantum field theory is described by the 
1-particle state vector in the Fock space 



(Pb|P'b) = (27r) 3 2u!B5 3 (P-P'), w B = J M 2 B + P 2 . (2.11) 



1 We choose the system of units where h = c = 1. Any 4-vector has the components a = (do, a), 
and the metric is = diag(l, —1, —1, —1) 
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However, there is no interpolating field in the Lagrangian corresponding to the bound state 
particle. The completeness condition of the Fock-state vectors in the presence of bound 
states reads 

r r/ 3 P 

1 = / |Pb >W< PbI + '" < 2 ' 12 ' 

where ellipses stand for the contributions of the states with elementary particles and from 
the multi-particle scattering states. 

Using the completeness condition (|2.12|) , it is straightforward to single out the bound- 
state contribution in the Green function ( |2.1| ) when P 2 —>■ M B (equivalently Pq — > w 2 B ). 
The quantity (p\G(P)\p') exhibits the pole behavior at this point 

ip\G(P)\p') = i + (P\R{P)W) , (2.13) 

where (p\R(P)\p') denotes the regular remainder of (p\G(P)\p') at the bound-state pole that 
emerges from the contribution of other states in the sum over Fock-space vectors. Further, 
(p|$P B ) stands for the BS wave function of the bound state 

(p|$ Pfl ) = $ P » = jdx e*" (0|T^(/3x)V 2 (-ax)|P B ) , 

($ PB b'> = $p>') = Jdx a*' (P B \TM(3x)M-ax)\0} = $ Pfl (p') <y<p ® 7 f (2.14) 
The bound-state equation that can be derived for the state vector |$p B ) by substituting 



Eq. ( 2.13 ) in the BS equation for the Green function (|2.2|) , formally resembles the nonrela 



tivistic Schrodinger equation for two fermions 

G-\P b )\^ Pb ) = 0, (V Pb \G-\P b ) = 0, with G \P)-G- 1 (P) = K(P), (2.15) 



or 



|$ Pb ) = G Q {P B )K{P B )\^ B ) , ($ Ps | = (<j> PB \K(P B )G (P B ) . (2.16) 

Here P B = (w B , P B ). Explicitly, in the momentum space, we arrive at the following equation 
for the bound-state wave function HI 



<S> PB (p) = Go(P;p) J 0^K(P;p,p')^ PB (p), 

*Pb(p) = J0ji ®r B (p')K(P;p', P )G (P;p). (2.17) 

This equation should be solved in order to obtain the mass M B of the bound state. It is 
obvious that both equations: for $ Ps (p), and for its conjugate $ Ps (p), lead to the same 
bound-state spectrum. 

Next, we derive the normalization condition for the BS wave function. To this end, it is 
useful to start from the following identity 
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G(P)G-\P)G(P) = G(P) => G(P)(G \P) - K(P))G(P) = G(P) 



(2.18) 



If P 2 is close to M|, one can neglect the contribution from R(P) in Eq. Q2.13|) . We substitute 
the latter into Eq. fl2.18|) , and perform the integration along the closed contour C that 
encircles only the bound-state pole at P Q = w B , in the complex P plane. 



c 



l$P S )($l 



c 



(Go\P)-K(P))dP 

(p + w B -toy{p -w B + io) 2 
dp 



l*p s )(^ 



(P Q + w B - iO)(P - w B + zO) 
From the Cauchy's theorem, one has 

f(z) dz 



(2.19) 



\n-l 



C (z-z s y 



±2m 



(n - 1)! dz 



^ - f{z) 



(2.20) 



z=z s 



where the function f(z) is analytic inside the contour C, and the choice of the ± sign 
depends on whether one integrates counterclockwise (+) or clockwise (— ) along the contour. 
With the use of the above formula, from Eq. ( |2.19| ) one readily obtains the normalization 
condition for the BS wave function 



d A p d A p' 



(2vr) 4 (2tt) 



d 

dPn 



(G^(P;p,p') - K(P;p,p')) 



Pq=w b 



2w B . (2.21) 



The equations (237] ), together with the normalization condition ( [2.21J) , completely determine 
the BS mass spectrum and the BS wave function. 

At the end of this section, we shall consider in some detail the spin content of the BS 
wave function. In particular, we shall demonstrate that one can rewrite this equation in 
terms of "fermion-antifermion" rather than "two fermion" wave function. 

We work with the following representation of Dirac 7-matrices 
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7o 







7 = 7 



a 



a 

(7 



-0:7 



(2.22) 



The free two-fermion Green function given by Eqs. ( |2.9| ), fl2.10Q , can be written as 



G (P;p) 



Gaa{Pl, 

G ba {pi) 



Gab{Pl) 

G bb (pi) 



G aa {P2) G ab (p 2 ) 
G ba {p2) G bb (p 2 ) 



(2.23) 



where G uv (pi), u,v — a,b is the 2x2 matrix (operator) in the spin space of the z-th particle. 
Further, the BS wave function of the two-fermion system can be written as a column 



$p s (p) 



{$aa(P;p)\ 

<5> ba (P;p) 
\$ bb (P;p)/ 



(2.24) 



where, again, the components $ U1J (P; p), u,v = a,b are the 2x2 matrices in the spin space 
of two fermions. 
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Now, it is straightforward to ensure that the BS equation (|2.17|) can be rewritten in 
terms of "fermion-antifermion" wave function \l/p B (p) 

= SM( Pl ) J -0- K(P;p,p')^ FB (p') S^(-p 2 ) . (2.25) 

The wave functions \l/p s (p) and $p s (p) are related by (see 0) 

*p B (?) = ( !H5 5 ! a ^. ; ^ c = -i( t ab & p } a :_ 5? ) > (2-26) 

where 



$ 6a (P;p) $ 66 (P;p)y V*bb(P;p)(X2 $6a(P;p)(X2 



C = W = {X T) < 2 - 27 > 

denotes the charge conjugation matrix. 



III. THREE-DIMENSIONAL REDUCTIONS OF THE BS EQUATION 

One of the reasons why the three-dimensional (3D) reduction of the BS equation is nec- 
essary, is the absence of the usual quantum-mechanical probability interpretation for the 
wave function $ Ps (p) due to the dependence of the latter on the 0-th component of the 
relative 4-momentum. Further, in the presence of the confining interactions, it is extremely 
difficult to construct a "reasonable" kernel K in four dimensions that describes these in- 
teractions - we are not aware of any, completely successful attempt. On the other hand, 
the concept of static (3D) confining kernels that corresponds to an intuitively clear picture 
of infinitely rising potentials in the coordinate space, has been extremely useful in many 
semi-phenomenological applications to study, e.g. the characteristics of heavy quarkonia, 
etc. 

For this reason, below we shall mainly consider the static BS kernels (i.e. the kernels 
which do not depend on the cm. momentum P and the 0-th components of the relative 
momenta po, p' Q ) 

K(P;p,p') - K st (p,p') = -iV(p,p') . (3.1) 

In this approximation, there are still different versions of the 3D equations for the bound- 
state wave function. Below, we shall consider these versions in detail. 

A. The Salpeter equation [6] 

In the approximation ( |3.1| ), from Eq. ( J2.17| ) it is straightforwardly obtained 

$ PB (p) = G (P;p) J |^ st (p,p')$P fl (p'), 

® pM = I S^b(pO^(^p)G (P;p), (3.2) 
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and 



$p s (p)= / ^33$P fl (p , )V(p' I p)Go(P;p) I (3.3) 

where 

^b(p) = /^*p*(p)> *P B (p)=/f^P B (p), G (P;p) = J^:G (P;p). (3.4) 
At the next step, we introduce the projection operators 

AP(Pi) = W *t W , Mft)=« W Pi + m i 7j , < = 1,2 J (3.5) 
with the properties 

E ^ = 1 , A^A^ } = K«>^ , ^(P,)AS ±} = • (3.6) 

With the use of the following identity 

ii + m = |(p i0 + ^)aS +) ( Pi ) + ( Pl0 - w; i )Ai~ ) (p i )}7o i) , (3.7) 
it is straightforward to obtain 

Go(P; P) = { p Air '""- P t n - &^k« a 7f 

I P - wi - w 2 + «0 P + + w 2 J 

= [P - fc x ( Pl ) - ^ 2 (p 2 ) ^(A^Pi, pa) - AS 2 - } ( Pl , pa)) 7 J 1} ® 7? } , (3.8) 

where Ag iaa) ( Pl , p 2 ) = A^(Pi) ® A^ 2) (p 2 ). 

Now the Salpeter equation Q3.3|) in the cm. frame (P# = 0) can be written as 



[p -Mpi)-^(p 2 )]Wp) = n(p) | ^ 7 ( 1 )® 7 ( 2 )y(p,p')$ Mfl ( P '), (3.9) 

where 



n(p) = (A^CPi,^) - Ak- } (p a ,P 2 )) = ^ + . (3.10) 

2^1 2w 2 

Introducing the "frequency components" of the wave function according to 

$P B (P) = E ^T 2) (P) , ^S a2) (p) = AS 2 ia2) (pi, P2 ) $p s ( P ) , (3.11) 

the Eq. (|3.9|) can be reduced to the following system of equations 



[M B T K + w 2 )p^\p) = ±A^\p, -p) 7 ? } ® 7f / ^ V(p, p') $ Mb (p') , (3.12) 

with additional conditions 

*£?(p) = 0, , WP) = $£ B +) (P) + ^(p) • (3.13) 



The normalization condition can be readily obtained from Eq. (|2.21|) by using the ap- 
proximation ( |3.1| ) for the kernel, the relation between 4D and 3D wave functions fl3.2| ), and 
the decomposition of the wave function Q3.11|) , ( |3.13| ) 

/ Wf { \^\v)\ 2 - \^m b \p)\ 2 } = 2Mb ■ (3.14) 
Note that the wave function $m b (p) can be represented in a form analogous to ( |2.24| ) 



*M fl (p) 



foft(p) 
$6a(p) 
V $bb(p) J 



(3.15) 



The constraints ( |3.13| ) can be considered as equations for the components $ a fe(p) and $ fefl (p). 
The solution of these equations gives 

®ab = (miii>2 + m 2 w i y 1 {wi(cr (2) p 2 )l> aa - w 2 (cr (1) pi)$ fef) } , 

$ 6a = (m x w 2 + m 2 w 1 )- 1 {w 2 (o- (1) pi)l> oa - tui(<r (2) p 2 )$»)} . (3.16) 

For the "frequency components" $^ lQ2 - ) = A^ 2 ia2 ^(p)l> :r .j / (p), x, y — a, 6, we obtain the 
following relations 

$^ ±} = ±(2(m lW2 + m^))- 1 {(wi ± m 1 )(w 2 ± m 2 )$ aa - (<r« Pl )(<r (2) p 2 )i> 66 } , 
= ±(2(miw 2 + mawj))" 1 {(o- (1) Pl )(cr (2) p 2 )<l aa - ( Wl T m 1 ){w 2 T m 2 )$ bb } , 
= ±(2{m x w 2 + mawO)- 1 {( Wl ± mi)(<r^p 2 )<£ aa - (u> 2 T m2)(<x (1) Pi)$»} , 

$£ ±} = ±(2(miw 2 + m^O)- 1 {(w 2 ± m 2 )(<T (1) Pi)$«« - ( Wl =f roi)(o-< 2 >p 2 )$»} . (3.17) 



The normalization condition ( 3.14 ) can be rewritten as 

J {ztt) 6 m\w 2 + m 2 wi 

At the end of this subsection, we shall consider the existence of the one-body limit in the 
Salpeter equation. From the physical point of view, it is clear that if the mass of one of the 
particles in the two-particle bound state tends to infinity, the equation for the wave function 
should reduce to Dirac equation for the light particle with a given interaction potential. Let 
us check this property for the Salpeter equation assuming, e.g. that the mass of the first 
particle tends to infinity. In this limit, 



S 



mi — > oo =>- W\ — > mi, 7q — ► 1, hi mi . (3.19) 

Then, the Salpeter equation for the bound-state vector \®m b ) is reduced to (E 2 = M B — 
mi) 

(E 2 -h 2 ) \® E2 ) = ± (l + ^) 7 f fe) . (3.20) 

Due to the presence of the prefactor | this equation differs from the Dirac equation 

for the particle 2 moving in the potential V - that is, the Salpeter equation does not possess 
the correct one-body limit. 

Now there arises an important problem to solve. We are willing to obtain the 3D reduc- 
tion of the BS equation in the static approximation, that correctly reproduces the dynamics 
of the system in the one-body limit - this property might be important for the description, 
e.g. the heavy-light qq bound states. 

Below, we shall consider several versions of the 3D reduction procedure, which lead to 
the correct one-body limit. 



B. The Gross equation [7] 

In the derivation of the Gross equation, first we assign a = and /3 — 1, in the definition 
of the cm. and relative momentum variables (|2.3| ). Physically, this means that the whole 
cm. momentum is carried by the particle 2. The free Green function has the form 

G (P;p) = -{^ + m 1 )®{P- p'+m 2 ) g (P;p) , 

9 ^ P] ^ = p 2 -m 2 i+ tO (P- P y-m 2 2 + t0 ' (3 ' 21) 
The first propagator can be rewritten as 

1 P = 1 2 - i^{p 2 - m\) , (3.22) 



p 2 — mf + iO p 2 — mi 

where the symbol P stands for the principal-value prescription. The approximation that 
leads to the Gross equation, consists in the substitution 

1 - -2«i 5(p °- Wl) . (3.23) 



p 2 — m{ + iO 2wi 

This approximation is called the "spectator approximation". Note that, in this approxi- 
mation it is not only the principal-value term in the propagator of the first particle that is 
neglected, but also the term containing S(po + w\) that emerges from 5{p 2 — m\). Conse- 
quently, in this approximation the particle 1 always stays on its mass shell defined by the 
equation p = w\. In a result of this approximation, the free Green function in the cm. 
frame (P^ = (Pq, 0)) can be rewritten in the following form 



9 



G GR (P ;p) = 2m5(p - w 1 )G^ ti (P ; p 



GR/ 



2iri8(p — Wi) 



Ai +) (P) 



Aj +) (- P ) f A'r'(p) : .\r'(-p) 



(+)< 



P - - u> 2 + zO P - Wi + w 2 + iO 



(3.24) 



where the functions G GR (Po;p) and G GR (Po; p) are related by Eq. 

After substituting Eq. Q3.24Q into ( 2.16|) and integrating over the variable po, in the cm. 
frame (now P M = (M B ,0)) we arrive at the Gross equation for the 3D bound-state wave 
function 



$M fl (p) = G GR (M B ;p) 



d 3 p' 

(2vr) 3 



V(p,p')$ MB (p') 



$M s (p) 



d a 



p 



(2*y 



®M B (p')V(p',p) Gq (Mb] p) • 



(3.25) 



Now, using again (|3.24 ) together with ( [3.6|) , we arrive at 



[M B - h{p) - h 2 {-p)\ $ Ms (p) 



(2tt) 3 



nP.P')*M fl (p'), (3-26) 



which has the correct one-body limit when mi — > oo. 

The normalization condition for the 3D wave function that satisfies the Gross equation, 
can not be obtained in a standard manner, by using Eq. ( 2.21 ). In order to demonstrate 
this, note that according to Eqs. (|2.16|) , (|3.24|) and ( |3.25|) , 4D and 3D wave functions are 
related by 

$ A / B (p) = 2n5(p - Wi) $m b (p) , ®mM = 2tt5(p - w x ) Im b (p) • (3.27) 

Now if in the normalization condition (|2.21 ) with the static kernel ( |3.1| ), the rela- 
tion ( 3.27 ) between the 4D and 3D wave functions is substituted, one arrives at the ill- 
defined expression containing the product of (^-functions with the same argument. For this 
reason, instead of the rigorous derivation, from the analogy with the Salpeter equation, one 
merely assumes that the solutions of the Gross equation satisfy the following normalization 
condition 



d 3 p 

(2tt) 3 



P 



I^(P) 



! } = 2M J 



B ■ 



(3.28) 



C. The Mandelzweig- Wallace equation [8] 

In the derivation of the Mandelzweig- Wallace (MW) equation, the parameters a and (3 in 
the expression of the cm. and relative momenta ( |2.3| ) are defined according to Wightmann 
and Garding 

a = aW = i±^g , = 0(s) = S - m l +m \ , = J * = l? _p». (3.29) 
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In the cm. frame, from Eqs. (|2.3| ) and (|3.29|) it follows 

/ ^ x / ^ ^ Ml, + m? — ml 
Pi = (Ex+p ,p), p 2 = (E 2 -p ,-p), E 1= m, E 2 



M B — m\ + m\ 



2M f 



2M, 
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E 1 + E 2 = M B , Ex- E 2 



2 2 

m\ — m 2 



M, 



B 



(3.30) 



Further, we define in the cm. frame 



dp 



G^(M B ;p) = / -^[Go( Pl ,p 2 )+G (pi,j9 c 2 r )], 



p^ = (E 2 +p Q ,-p), (3.31) 

where Go(pi,p 2 ) is given by Eqs. Q2.8|) , (|2.9| ) and (|2.10|) . After integrating over p , we obtain 

Ar2 + (P,-P) A+-(p,-p) 



E x + E 2 - wi - w 2 + iO + -Ex + E 2 + wx + w 2 



Ai2 + (P, -P) 



A12 (p, -P) 



+ 5 7 W (g, 7 2 ) . (3.32) 

Ei - E 2 + Wl + w 2 Ex + E 2 + wx + w 2 ) !0 !0 v ; 

The MW equation is obtained from the BS equation in the static approximation, by using 
the combination G {px,p 2 ) + G (px,p 2 ) instead of G (px,p 2 ) alone. Unlike the Salpeter 
version, now all four possible projection operators A^ , A^ - '', A-rz + \ ^-V2~\ enter the 
expression of Gfl™ (M B ;p), Eq. ( |3.32| ). For this reason, the inverse operator for the free 
Green function in the 3D space exists, further, in analogy with Eq. (|2.15 ), we can define 
the inverse of the full Green function in 3D space according to 

[GrV\M B - p, p') - [G MW }-\M B ; p, p') = V(p, p') , 



where 

G MW (M B ; p, p') = (2vr) 3 5 3 (p - p') Gtf"(M B ; p) . 
The MW equation for the bound-state vector \&m b ) * s given by 



(3.33) 
(3.34) 



[G 



MWl-l 



<W = o 



\<s>M B ) = Grv\$ MB ) 



(3.35) 



Note that we can rewrite the inverse of the free Green function in the MW equation as 



[G 



rr 1 =TS 1, ®7o ca) 



(Ex -h x )® — + {E 2 -h 2 )® — 

W 2 Wx 



With the use of this identity, one can rewrite the MW equation as 

\*m b )='yP®'y!? ) V\$ Mb ). 



(Ex -hx)® — + (E 2 -h 2 )® — 

W 2 Wx 



(3.36) 



(3.37) 



Let us now consider the limit of this equation when mi — > oo (see Eq. (RTT^)). In this 
limit, according to Eq. ( |3.30|) , Ex — > mi, E 2 — > M B — mi, and the equation ( |3.37|) simplifies 
to the Dirac equation 



[E 2 -h 2 ]\§ MB )= 1 f ) V\§ MB ). 
Consequently, the MW equation has the correct one-body limit. 



(3.38) 
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D. The Cooper- Jennings equation [9] 

The parameters a(s) and (3(s) in the Cooper- Jennings (CJ) version are chosen as 

a(s) = , ° l(a) , , , fts) a2(s) 



oti(s) + a 2 (s) a 1 (s) + a 2 (s) 



s + mj-m 2 2 s-ml + mt, 
ai ^ S ' = Vs ' = 2~Js ' ^ ' 

The free Green function for the CJ equation is given by 

G^(P;p) = -(jfi + m a ) ® (p' 2 + m 2 )g^(P;p) , (3.40) 

where gQ 3 {P',p) is constrained by the elastic unitarity and can be written in the following 
form 

gZ J (P;p) = 2m f° ^li^l 5+[(a(s')P' + pf - m\\ 5+[((3(s')P' - p) 2 - mg] . 

J{m 1 +m 2 ) 2 S ' — S — l\J 

(3.41) 



Here 5 + (x 2 — a 2 ) = (2a)~ 1 5(x — a), P' = ys'/sP, and the function f(s,s') satisfies the 
condition f(s,s) = 1. 

After integration, the expression (|3.41|) yields 

g™(P;p) = -2J-^- ^ss~ p l^A ai (s p ) a 2 (s p ) , (3.42) 



where s p = (\Jm 2 — p 2 + yrnj — p 2 ) 2 . 

Choosing the function f(s, s p ) in the form 



ss p — [m\ — m 2) 

we arrive at the following expression for g^{P\p) 

PT , N 2s 2J~s5(2P-p) 
g^(P-p) = -2m - - - — i ^ . 3.44 

Note, that in the cm. frame, 2y/s6(2P ■ p) = 5(po). Because of the presence of the 5 
function, one can rewrite the free Green function from ( |3.40| ) in the following form (again, 
in the cm. frame) 

G^(M B ;p) = -(^ + mi )®^ 2 + m 2 )g^(M B ;p) } p 1 = (E 1 ,p), p 2 =(E 2 ,-p), (3.45) 

where E\ and E 2 are given by Eq. ( |3.30|) . The free Green function for the CJ equation in 
3D space is related to 4D Green function according to 

G^(M B ;p) =27n5(p )G^(M B ;p), (3.46) 
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where Gq j (M b ; p) is given by 
G CJ (M B ;p) 



, ;J/M ^^ _ 1 ($1 + mQ g (#2 + ^2) _ 1 (#1 + g ($2 + m 2 ) 
2{w 1 + w 2 ) Ej-wj ~ 2(w 1 +w 2 ) E\-w\ 



1 ^ + 777,1 1 $ 2 + m 2 



2(wi + w 2 ) #2 - ra 2 2(w\ + w 2 ) $1 - mi' 
In the limit , when one of the masses tends to infinity, 



(3.47) 



2(wi + w 2 ) 



at rrii — > 00, i = 1,2. (3.48) 



Consequently, the CJ equation has the correct one-body limit. 

Note that, using the properties of the projection operators, the free Green function in 
the 3D space can be rewritten in the following form 

G^(M B ; p) = ] {( Wl + EJfa + E 2 )A ( + + \p, -p) 

£\W\ -+- W 2 )CL 

- (wi + E t )(w 2 - E 2 )A { +-\p, -p) - ( Wl - E 1 )(w 2 + E 2 )a[ 2 + \p, -p) 
+ (wi - E 1 )(w 2 - E 2 )A[ 2 -\p, -p)] 7^ ® 7? } , (3.49) 

where 

a = El-wl = El-wl = -\Ml + bl-2{wl + wl)}, b = E x - E 2 = . (3.50) 

4 IVl r 



E. The Maung-Norbury-Kahana equation [10,11] 

The free Green function for the Maung-Norbury-Kahana (MNK) equation is again given 
by Eq. flOCD , but with 

5 + {[(a( S )P + p) 2 - ml} 1 -^ - l(P(s)P - pf - ml) 1 -^} 
9o (P; P) = -2« [( a ( s )p + pj2 _ m 2] + [ (/3 ( s )p _ p y _ m 2] + 20 > ( 3 - 51 ) 

with 7/ = (mx ~~ m 2 )/(mi + 777.2). This Green functions, of course, satisfies the unitariry 
condition in the elastic channel. In addition, it has the property that the particles 1 and 
2 in the intermediate states are now allowed to go off mass shell inverse proportionally to 
their masses - so that, if one of the particles becomes infinitely massive, it is automatically 
kept on its mass shell. 

After some transformations, the Green function from Eq. (|3.51 ) in the cm. frame can 
be rewritten as 

MNK, w x _ 2 • $(Po - Po) + _ R-b _ + , 



i? = Jfe 2 - 4?/ 2 a , 6 = M B + b y , (3.52) 
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and 60 is given by Eq. ( |3.50|) . By using the above expression, we obtain 

G (Af B ,p) = 9B / + - , \ > Pi = (^i+Po.P). = (£2 -Po,-p). (3.53) 

zrt^p p -+- a) 

The relation between the 4D and 3D free Green functions in the MNK version is given 

by 

G^ K (M B ;p) = 2iti8(p Q - p+) Gf NK (M B , p) . (3.54) 

Using the properties of the projection operators, the free Green function of the MNK equa- 
tion can be recast in the following form 

GT K (Mb, P) = 9g , Jp- , {[(wx + E 1 )(w 2 + E 2 )A[ + 2 +) (p, -p) 
ZK(p p + a) L 

- ( Wl + Ei)(w 2 - E 2 )A[t\p, -p) - (w 1 - E 1 )(w 2 + E 2 )A { 12 +) (p, -p) 
+ ( Wl - E 1 )(w 2 - E 2 )A { h -\p, -p)] 

- [p+po + (wi - ^2)p^(a£ +) (p, -p) - Air } (p, -p)) 

+ K + ^(A^p, -P) " Air } (P, -P))]}7^ ® 7^ • (3-55) 
In the limit when mi —* 00, the function gQ SNK (P',p) from Eq. ( |3.51|) is reduced to 

So MNK (^P) 



— 2?ri 5(po) 

an — x p\ TTl 2 lTTl\ 



2 „,2 o... ' ( 3 - 56 ) 



From Eq. ( |3.40 ) we can evaluate G^ NK (P;p) in this limit: 



- 2m *M + + j (3 . 57) 

mi— >-oc 2uii p 2 17l 2 



where p\ and p 2 are defined by Eq. ( |3.45| ). Integrating this relation over p , for the 3D free 



Green function in the cm. frame we obtain 

mj^oo 27711 p 2 TTl 2 



G (M B ;p y 



,2 ,,,2 • (3-58) 



Since the factor (pi + mi)/(2mi) tends to unity in the limit m\ — > 00, one concludes that 
the MNK equation has the correct one-body limit. 

F. The normalization condition for the wave function in MW, CJ and MNK versions 

The 3D free Green function in either of MW, CJ, or MNK versions, in the cm. frame 
can be rewritten in terms of the projection operators: 

CoW,.P) = E ± g '7S ;P> Alr 2> (P. -P) 7^" ■ P=|P|, (3.59) 
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where 



MW : D {aia2) = - ^ — — , d = 1 , 

+ w 2 ) - [oLiEx + a 2 E 2 ) 

CJ : D {aia2) = (E 1 + a 1 w 1 )(E 2 + a 2 w 2 ), d = 2{w l + w 2 )a . 
MNK : £>(«!<») = (Si + a l w 1 )(E 2 + a 2 w 2 ) 

R-b (R-b 



2y V 2y 



+ (Ei + aiiui) - (S 2 + a 2 u>2)^ 



d = 2RB, B = -(^— - + b ) +a, (3.60) 
2y V 2y J 

with Si, S 2 , a, 6 , -R5 2/ defined above. 

The equation for the bound-state wave function frequency components ( |3.11| ) can be 
directly obtained from Eq. ( p.3[ ) by substituting the above expression for the free 3D Green 
function and using the properties of the projection operators 



[M B -(a lWl + a2^)]$ir 2) (P) =^ (aia2) (M B ;p)ASr 2) (P,-p)7S 1) ®7S 2) x 

(2tt 



x E /^^(P.P0^V), (3-61) 



*1"2 



where, for the different versions 
MW : A (±±) = 1 , A (±T) 



M B 



Wi + w 2 



C J . A {a W ) = M B + (<*lWl + U2W 2 ) 

2(wi + w 2 ) 



MNK : A^ a2) = — *— I a [Mb + (aittfi + a 2 ^ 2 )] - [M B - (a^i + a 2 w 2 )]^—^- x 

ARB L 



2y 



x 



+ (Si + a x w x ) - (S 2 + a 2 w 2 M . (3.62) 



Thus, the MW, CJ and MNK equations couple all four frequency components of the wave 
function: &Mg \ > ^Mb anc ^ ®Mb ■ ® ne can f° rma lly extend these notations for the 

Salpeter (SAL) and Gross (GR) versions, defining 

SAL : A {±±) = ±1 , A i±T) = , 

GR : A i+±) = ±1, A { - T) = 0. (3.63) 

It is immediately seen that Salpeter and Gross equations couple only two frequency compo- 
nents of the wave function, other two being equal to 0. 

For the derivation of the wave function normalization condition in MW, CJ and MNK 
versions, let us consider the full 3D Green function that obeys the equation 



15 



& = G l + GiyG 1 = & + &V& Q , % = MW, CJ, MNK . (3.64) 
In analogy with Eq. ([2.131) , this Green function develops a bound-state pole(s) 



(p|G(P)|p') = £ 



(p|$P fl )($ Ps |p / ) 



b P 2 ~Ml 



(pI^(p)Ip'). 



(3.65) 



This leads to the normalization condition in MW, CJ and MNK versions in analogy with 
Eq- (TO 



d 3 p d 3 p' =< 



<9 



<9M, 



B 



(G^Msjp,?'))- 1 - y(p,p') $V,(p') = 2M B . (3.66) 



(2tt) 3 (2tt) ; 
Now using the relation 

(& (M B ; p, p'))- 1 = (2n) 3 5 3 (p - p') (6* (M fl ; p))" 1 



(3.67) 



and the fact that V(jp,p') does not depend on M B , the normalization condition can be 
rewritten as 



I $ M B (P) 



(2vr 



5 



(G^Mbip))" 1 $ j Mb ( P ) = 2M b . 



(3.68) 



If now one substitutes here the expression for the free Green function given by Eq. ( p.59 ), 
one obtains (below, we drop the superscript "i" labeling various versions) 



CtlCt2 



(2tt 



(3.69) 



where 



ft a2) (M B ; P ) 



d ( d(M B ;p) 



dM B \D a ^(M B ;p) 

By using the explicit expressions for £)( aia2 ) and d given by Eq. fl3.60| ), we obtain 

?(«i«2) _ + a 2 E 2 



(3.70) 



MW 
CJ 
MNK 



12 



M. 



B 



f {ai.cv) _ 2 ( w i + w 2) a x WxE l + a 2 w 2 E 2 

J 12 — 



M B (Ei + aiWi)(E 2 + a 2 w 2 ) 



AaiOL2) 
J12 



B 

]j{aia2) 



M B B , 9x Ml 



R-M B \ 2 



M b + 



aiWi + a 2 W2 



i? 



M 



B 



+ 2 



2?7 

R-M B ^ 2 

2y 



U 'R-M B M B y\ 
(a 1 w 1 - a 2 w 2 ) I — h 



2y 



J 



(3.71) 
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Let us note that the normalization condition ( |3.69|) is valid for the Salpeter and Gross 
versions as well, provided we choose 

SAL : A ( r 2) = ^^, 

GR : A ( 2 +±) = l, /£- ±} = 0. (3.72) 

Let us emphasize that the Salpeter, MW, CJ and MNK equations can be used for the 
bound systems with the equal masses of the constituents, whereas the Gross equation can 
not - the particle "1" (spectator) should be heavier than the particle "2". This is due to 
the approximation (|3.23|) that was done in the free Green function of the Gross equation. 



Further, it is directly seen from Eqs. (|3.62| ) and ( |3.63| ), that the Salpeter and Gross equations 
are linear eigenvalue equations for determining Mb (the functions A^ ai ° 2 ^ do not depend on 
M B ), whereas MW, CJ and MNK equations are not, and Mb enters the right-hand side of 
these equations as well. 

Let us now concentrate on the properties of the coefficient functions in detail. In the 
case of the equal-mass constituents mi = m 2 = m and W\ = w<i = w, we obtain 



MW : A {±±) = 1 , A (±T) 



Mb 
2w 



CJ : A<**> ~ MB ± 2W , A^) = Ml 

Aw Aw 



MNK : A^ = M b±Z™ A ^ = -. (3.73) 

2Mb 2 v ; 

It is immediately seen that in the equal- mass case, the Mb drops out from the equations 
for mixed components $mJ^ in the MW and CJ versions - that is, these components are 
redundant and can be eliminated in this case. 

The functions /i 2 1Q2 ^ in the equal- mass case are given by 

MW : ^ = ±1, A ( 2 ±T) = 



(M B ± 2w) 2 ' 



CJ : fu — ^JTJ r~o \o > fu 



MNK A±±)_ 2{{M B ±2 W f-^) (±T) _ 

MINK . / 12 - {Mb±2w)2 > /« ~ 2 - ^- 7A ) 

From these expressions, we immediately see that in the CJ and MNK versions the func- 
tion ^ has the second-order pole at 



M B -2w(p s ) = 0, Ps = - yjM B - Am 2 . (3.75) 

It can be shown that in the non-equal mass case in the CJ and MNK versions the function 
f[ 2 ^ has the second-order pole at 



Ef - w 2 ( Ps ) = , Vs = \^Ml+ (^^) 2 - 2K + m 2 ) 2 , (3.76) 
while the other components f^ + \ fi2 and f-fe do not have any poles. 
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G. Logunov-Tavkhelidze quasipotential approach [12] 



There exists the theoretical possibility to construct the 3D analogue of the BS equation 
without using the instantaneous approximation. To this end, one may use the Logunov- 
Tavkhelidze quasipotential approach formulated in Ref. \12\ for the case of two spinless 
particles, and generalized in Ref. Jl3| to the case of two fermions. 

We introduce the following definition. For any operator A(P) in the momentum space, 

(plI(P)lp') = / (P\MP)W) • (3-77) 

Then, from Eq. ( |2.2| ) one obtains 

G = G + G^KG , (3.78) 

where Go is given by Eq. (|3.8|) . 

Due to the fact that the operator II defined by Eq. ( |3.10| ) can not be inverted, the inverse 
operator of Go does not exist as well. As a result, one can not define the interaction potential 
by the formula analogous to Eq. (|3.33| ). In order to overcome this problem, it is convenient 
to introduce the Green function Go defined by 

Go(P; p) =i[P Q - MPi) " ^ferSo 1 } ® 7? } , 
G (P; P , p') = (2tt) 3 5 3 ( p - p') Go(P; p) 7 « ® 7 <) 2) ft (3.79) 

where ft = Hjo^ <8> ■ 

Now, the inverse of the operator Gp_(P; p, p') = (27r) 3 <5 3 (p — p') Go(P; p) exists, and one 
may define 

G = Go + Gol<G, (3.80) 

from which follows that 

G = G + Go(l-7S 1) ® 7 S 2) ft). (3.81) 

It is clear, that near the bound-state pole the Green functions G and G differ only by 
the regular term, since Go is regular in the vicinity of the pole. Consequently, in order to 
derive the bound-state equation, one may use G instead of G, and define the interaction 
potential according to 

[Gor'-iar^v, (3.82) 

from which it follows that 

G = Gp_+GoV:G, (3.83) 

and 

V=[G ]- l G Q l<G\G G }- 1 . (3.84) 
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For any given kernel K(P;p,p') the interaction potential can be constructed by using 
Eq. ( p. 84 ). The equation for the bound-state wave function in the cm. frame can be 
obtained directly from Eq. ( |3.82| ) 



G \M B ) 



\®Mb) =°> 



1$ 



GV\^ Mb ) 



Defining the quasipotential as 

V g (M B ;p,p') 

we obtain 



■ (i) 
^7o 



M B - hi(p) - h 2 (-p) 



M b ) 

4 2) V(M B ;p,p') 

d 3 p' 



(2k) 



V g (M B ;p,p') $ MB (p') 



(3.85) 
(3.86) 

(3.87) 



Note that in the instantaneous approximation ( |3.1| ), the interaction potential reduces to 

V(M B ; p, p') = 7 S 1} ® 7 f n(p) 7 S X) ® 7 f K at (p, P') • (3-88) 

As a result, the quasipotential equation reduces to the Salpeter equation. 

The first-order quasipotential is defined by Eqs. ( |3.82|) and ( ft.86| ), if in the former the 
full Green function G is substituted by the free Green function Go 

l/«(M B ;p,p') = z 7 ( 1) ® 7 f (pllGvr'G^KGlGo]- 1 ^) (3.89) 

from which, in the static approximation one obtains 

VU(M B ; p, p') = n(p) 7 « ® 7 f z*T st (p, p') . (3.90) 

It is seen that, unlike the full quasipotential equation, the first-order equation does 
not reduce to the Salpeter equation in the static limit. Only when one may neglect the 
negative-frequency component of the bound-state wave function, the first-order equation 
again reduces to the Salpeter equation in the static limit. Here we note, that the first-order 
quasipotential equation was used in Ref. |14| in order to evaluate the dynamical retardation 
effect in the qq bound system mass spectrum (i.e. the effect that stems from the deviation 
of the BS kernel from the static one). 

In the rest of this subsection, we consider the normalization condition for the quasipo- 
tential bound-state wave function. Near the bound-state pole, the 3D Green function G(P) 
develops a pole (|3.65|) . Using the fact that in the vicinity of the bound-state pole the Green 
functions G(P) and G(P) coincide up to the regular term, it is straightforward to obtain 
the normalization condition 



d 3 p d 3 



P' 



(2k) 3 (2tt) : 



$M fl (p) 







dM, 



B 



(Go(M B ; p, p'))- 1 - V(M B ; p, p') ) $ Mfl (p') = 2M B . (3.91) 



From this equation, using the definition of the conjugate wave function (|2.14j) and Eqs. ( |3.79| ), 
( |3.8fi| ), we obtain 



d 3 p d 3 p' ~ 



(2vr) 3 (2k) 3 



$L(P) 



1 - 



d 



dM 



B 



V q (M B] p,p'] 



M B (P 



2M, 



B ■ 



(3.92) 



As it is seen from Eq. fl3.92p , in the static limit the above normalization condition 
reduces to the normalization condition for the Salpeter wave function only if one neglects 
the contribution from the negative-energy component of the wave function. 
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IV. MESON SPECTROSCOPY 



A. Partial-wave decomposition 



The properties of the qq bound systems in the 3D formalism obtained from the BS 
equation in the static approximation, were studied in Refs. [TT|JT^-^5|, without making 
any additional assumptions. Note that these 3D equations can be written either, as in 
Eq. d^), for the 2-fermion bound-state wave function [|g|T3]^J2T|^^,0,§3]§§ , or for 



the fermion-antifermion bound-state wave function p0| , p4| , p5| , p8| - |32]j3^| (the latter is obtained 
from Eq. ( |2.25|) in the static approximation.). Further, one may write down these equa- 
tions it terms of either the frequency components of the 3D wave functions 1^' (p) and 



if?MB^ (p) (^ ne latter denotes the frequency components of the fermion-antifermion wave func- 
tion), or on terms of their linear combinations $ aa (p), $ feb (p), etc, \I/ aa (p), ^^(p), etc, see 
Eqs. ( |2.24j ), ( |2.26| ). Below, we shall use the form of the 3D equations given by ( ^.61[ )-( pi33] ), 
with the normalization condition given by ( |3.69| )-( |3.72D flIT| , |2"6| , |3"5 



In order to rewrite the equations explicitly in either of the forms above, one has to specify 
the explicit spin structure of the interaction potential. This potential consists of several 
parts. First, there is the one-gluon (OG) exchange piece dominating at short distances. In 
the Feynman gauge, the spin structure of this piece is given by 7^ 



(i) 



7 



(2)ft 



7$ 1} ®7S 2) 



7W <g> 7^). In accordance with the static approximation, however, we neglect the second 
term in this expression 0. In addition, there is the confinement (C) piece in the potential 
that dominates at large distances and leads to the formation of the qq bound states. The 
spin structure of this piece is not known a priori. We choose it to be the mixture of a 
scalar and a zeroth component of a vector. Further, sometimes an additional "instanton- 
induced" piece corresponding to the t'Hooft interaction, is included in the potential |25| . 
The spin structure of this term is given by the equal mixture of scalar and pseudoscalar 
parts. The rationale for including the latter piece is the following. In the absence of the 
proper treatment of the Goldstone nature of light pseudoscalar bosons that is due to the 
spontaneous breaking of chiral symmetry in QCD, the t'Hooft interaction mimics this effect, 
leading to the large mass splitting between the pseudoscalar and vector mesons. Note that 
the chiral symmetry can be consistently incorporated in the 3D framework (see, e.g. [f2~0|j), 
albeit at a cost of the more involved formalism. For example, in this case the Hamiltonian 
of the free quark is replaced by 



hi(pi) = at^pi + ma^ 



S i (pi)a»p i + A(Pi)7o 



(0 



(4.1) 



where Aj(pj) and -Bj(pi) are determined by solving the gap equation for the quark propagator 
with the static potential. Below, however, we do not consider this approach. 
Thus, the spin structure of the static potential we shall be using, is given by 



V = 7 « ® 7 J 2) Voc( r ) + (^7o 1} ® l { o ] + (1 - ® I {2) ) V c (r) 

+ (/(i) (g) /(2) +7 (i) (g)7 (2) )V T {r), 



(4.2) 



where the last term corresponds to the t'Hooft interaction, all potentials are assumed to be 
local, and < x < 1. 



20 



Let us now turn to the wave function. It is possible to "solve" the constraints imposed 
on the frequency components, defining 



toi+oimi 



/* Q2) (P)\ 
$i? Q2) (p) 

^r 2) ( P ) 

V^ 1Q2) (p)7 



(4.3) 



where 



^ ( 2 aia2) (p) 



Iwi + atimi w 2 + a 2 m 2 



2wi 



2w 2 



(4.4) 



and Xm 1&2 ^(p) * s the unconstrained Pauli 2x2 spinor. For this spinor, the following system 
of equations is obtained 

[M B - (a lWl + a 2 w 2 )] X ^ B a2 \p) = A^(M B ;p) £ / ^V^ a ^(p,p') X ^(p% 

a' ia ' 2 J ^1 

(4.5) 



where 



+V T (p - p') (B}' 
and 



Klr^^P, P) = ^ ia2) (p)(^i(p - p')5i aia2CW (p, P') + V 2 (s; P - p')Bp a ^(p, pO 



(P, P) - ^° 



(p,p')))^ Q ^(p') (4.6) 



51 



(aiaja'jaj) 



P,p') 



B 



(aia2«' 1 «2) 



P,P 



ll, p)(ff l2, p)(ff ( ¥)(ff (J) pO 

(wi + a 1 m 1 )(w 2 + a 2 m 2 )(w[ + a[mi)(w2 + a 2 m 2 ) 
aia[ (er 

WpX^V) a 2 a' 2 (^ (2) P)^ (2) P') 



(wi + aim 1 )(w' 1 + a'imi) (w 2 + a 2 m 2 ){w' 2 + a' 2 m 2 ) 
«p)(<r( 2 )p) ^(o-Wp'X^V) 



(wi + aim 1 )(w 2 + a 2 m 2 ) (w[ + a! 1 mi){w 2 + a' 2 m 2 ) 

OL\0?2 ( CT 

( 1 )p)(cr( 2 )p') 

(wi + aimi)(w 2 + a 2 m 2 ) (wi + a; / 1 mi)(iy 2 + a 2 m 2 ) 



(4.7) 



Vi = K) G + Vc, V 2 (x) = Vb G + (2x-l)V c . (4.8) 

The functions Vog(p — p')) ^c(p — p') an d Vr(p — p') are the Fourier-transform of the local 
potentials Vbc( r ); ^c( r ) an d Vr(r), respectively. 

The normalization condition for the Pauli spinors Xms° 2 ^(p) follows from (|3.69|) 

E / 7^ Xlt ia2 \p) ft a2 \p) xlr 2) (p) = 2Mb • (4.9) 



21 



The partial- wave expansion of the Pauli spinor Xm^ (p) * s given by 

xfe^PH E xtri(p)= E (nlLSJMjB&^ip), n=£, (4.10) 

LSJMj LSJMj P 

where Rlsjmj(p) denote the radial wave functions, and L, S, J, Mj stand for the total orbital 
angular momentum, total spin, total angular momentum, and the projection of the total 
angular momentum of the qq system, respectively. 

The partial-wave expansion of the potentials reads as 

V(p - p') = (2vr) 3 ]T (n\LSJMj)V; L (p,p')(LSJMj\v!) , I = OG,C,T, (4.11) 

LSJMj 

where 

2 r°° 

Vi L (p,p') = - r 2 drj L (pr)Vi(r)j L (p'r), (4.12) 

7T JO 

jl being the spherical Bessel function. 

Using the fact that for the spherical potentials Vi(p — p') = Vj(\p — p'|), one may write 

V iM) = 4^2 j\ dzP L(z) V l ^p2+p> 2 -2pp>z) , (4.13) 

where Pi(z) denotes the Legendre polynomial. The above form is convenient when the 
function Vi(p,p ; ;z) can be written in the analytic form. 

In order to carry out the partial-wave expansion in the bound-state equation, it is con- 
venient to introduce the operators S = \(<J^ + cr^), cr = \{<r^ — cr^), instead of the 
individual spin operators o-M, i — 1,2. At the next step, one uses the known values of 
matrix elements of the operators Sn, crn, and the tensor operators 

S 12 = 3( CT «n)( CT ( 2 M - (<r«er( 2 )) = 6(Sn) 2 - 2S 2 , (a^a^) = 2S 2 - 2 (4.14) 

between the different spin-angular momentum states 

(LSJMj\ (j^) \L'S'J'Mj,) = 6jj,5 MjMjl (L\\n\\L')(S\\ ) \\S') x 

xW(LL'SS';U)(-lf +L - J , 
(L||n||L') = V2L' + 1 (L'100|L0), 

(S\\S\\S') = 5 SS ^S(S+1)(2S+1), 

(S\\tr\\S') = (-l) s+1 V36 S \s'-i\ 

(LSJMj\S 12 \L'S'J'Mj,) = 5jj,5 MjMIj 8 ss > { ^ LSJ \\^\^' J A = (_i)^ 51 x 

x ^120(2L' + 1) < L'200|L0 > W(LJ21; 1L')| , (4.15) 
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where W stand for the conventional Racah coefficients. The bound-state equations for the 
radial wave functions R^lsT\p) are then obtained straightforwardly 

[M B - (a lWl + a 2 w 2 ) ] ifc 2 >) = A^(M B ;p) ]T /°° p' 2 dp' x 
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+ (a 1 a[Xi 2 - a ^\p)Mr' ia ' 2 \p') + a 2 a'Xr- a2 \p)M2' 1 ' a2 \pW}J(^P,P')) ^>') 

- {a 1 a[N[ 2 a ^\p)U^\p') - a 2 a'X2 1 ' a2 \p)M2' 1 ~ a ' 2 \p'W2eAx ] p,p')) R^ip') 

(Mr 2 \pW}2 ia ' 2 \p) + a 1 a 2 aW 2 M2- ai - a2 \p)M2~ a ' 1 - a ' 2 \ P > )) 

± {a ia2 N[ 2 a ^\p) + aW^'^ip'))) V T J (p,p') 

- {a l a\U[ 2 a ^\p)U[ 2 a '^\ 

± (a.a^^ip) + aWXr'^ip'))) V^jM)) R^P') 
+ ((a ia [Mr ia2 \pWi2~ a ' ia2 \p) - a 2 a ! 2 Mt~ a2 \pW^~ a ' 2 \p) 



T {a^N^-^ip) - a'^Mr'^ip'))) V TeJ (p,p')) R^ip')] } 
[M B - {a lWl + a 2 w 2 )]R^j(p) = A^(M B ;p) £ f°° p' 2 dp' x 



(4.16) 
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(4.17) 



where 

vB(p,p) 

Vkqj{PiP') 
V A (j±i)j(p,p) 



2J+1 



J 



J 



J(J + 1) r 



2J+ 1 
1 

(2J+1) 2 



^A _1 (P)P') ~ Va +1 (.PiP') 
Vi ±l { P ^)+AJ{J+l)Vl ±l ){ P ^) 



A = 1,2,T. (4.18) 



The normalization condition in terms of radial wave functions has a particularly simple 
form 
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OlQ2 



P 2 dp ( ai aa) 



/^ lQ2; (M B ;p) 



1? 



(aia 2 ) 



(P) 
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2M B . 



(4.19) 



B. Dynamical input 

For solving the bound-state equation, one needs further to specify the inter-quark poten- 
tials Vog, Vc, Vt, introduced above. Let us start from the confining part of the potential. 
It is believed that the explicit form of this potential (i.e. its dependence on the inter-quark 
distance) is in principle, derivable from QCD. At present, however, the only tangible theoret- 
ical constraint on the form of this potential is the linear growth at large distances obtained 
within the quenched lattice QCD . Less compelling arguments based on the background 
field technique, were provided to justify the harmonic oscillator-type (~ r 2 ) behavior of the 
confining potential at small distances. With no rigorous solution of the problem in sight, 
one may use the potential that interpolates between the "known" behavior of the potential 
in different limiting situations P8|,|3!| (for a slightly modified version, see 



Voir) 



4 , 2 \ ( ^2^y 
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(4.20) 
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m l2 = m 1 +m 2 



fJ-12 



m 12 



(4.21) 



where Q 2 is the momentum transfer squared, and the factor | comes from the color- 
dependent part of the qq interaction, rif is the number of flavors (rif = 3 for u, d, s quarks, 
rif = 4 for u, d, s, c quarks, rif = 5 for u, d, s, c, b quarks). u , V , A , Aq C d are considered 
to be the free parameters of the model. The potential given by Eq. ( f4.20|) , effectively reduces 
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to the harmonic oscillator potential for the light quarks u, d, s, and to the linear potential 
for the heavy b, c quarks, that meets our expectations. In these limiting cases, the potential 
takes the form 



LINEAR 



HARMONIC 



Voir) 



r s(mi2) U 12 



\m\ra 2 



r — V Q I = ax + b x r 



V c (r) = ~a s (m? 2 ) r*-V )=a 2 + b 2 r* 



The one-gluon exchange potential is given by the standard expression 



Vbc r) = -- 



4 «s(m 



12 J 



b-ir 
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Noting that 



hm(-)"|- ' 
v^o dr] 



one can rewrite the potentials in the momentum space 

(9 2 / 47T 

LINEAR : Vc(p - p') = a 1 (2 7 r) 3 5 3 (p - p') + h lim — 

v-*o or] z V | p — P | 



HARMONIC 
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(4.22) 

(4.23) 
(4.24) 

(4.25) 
(4.26) 
(4.27) 



In order to investigate the properties of the qq bound systems, the linear potential was 



used both in the configuration space O,p0,[n],B3[], and in the momentum space PT| , |2"2"| , |3"3" ,4Ti 
In the latter case, a special numerical algorithm based on the the regularization ( |4.25|) , was 
utilized [ 2"2" . 4"0 ]. In Refs. P5] , PS| - |31]|J3^] , the matrix elements of Vc(r) were calculated in 
the configuration-space basis, in order to encompass the difficulties related to the singular 
character of the linear potential in the momentum space. 

The investigation of the qq systems in the framework of Salpeter equation was carried 
out |T7|- p0l , p6| , using the harmonic confining potential. MW, CJ and MNK equations with 
the harmonic confinement were considered in Refs. fTT| , |3"5|| . 

Some mathematical problems arise if the one-gluon exchange potential with the fixed 
coupling constant b-i is used for the calculation of the characteristics of qq bound systems. 
Namely, as it was shown in Ref. |p8,[30|, in this case the Salpeter wave function is divergent 
at r — > 0. For the running coupling constant this divergence is less pronounced but still 
present - now, the problem occurs in the decay observables which depend on the value of 
the wave function at r — > 0. In order to cure this divergence, in Refs. |]28|j30| the following 
regularization was proposed 



V 0G (r) = - 



4a s (r) 



V OG (r) = a a r + b, 



9 l 



for r > tq 
for r < r 



(4.28) 
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where 



2 ln(e-(^ a ) I a + e A ^ 2a ^] 



1-B 



ln(2 lm^T^/a + e 1/2 )) 



21n(e-' ia /a + e s / 2 ) 



(4.29) 



where a = Kqcd^-, 7 = 0.577215 • • • is the Euler-Mascheroni constant, and a sa t = 0.4, fi = 4, 
/i = 20. Further, 
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6(53 - 19n/) 
(33 - 2n/)2 



(4.30) 



Note, that in Ref . , the choice 5 = is adopted. The constants a g and b g from Eq. (|4.28|) 
are determined from matching of the potential and its first derivative at r = r . It turns out 
that the dependence of the qq system mass spectrum on the regularization parameter r is 
very weak provided the latter is chosen to be sufficiently small. 



C. t'Hooft interaction 



The t'Hooft interaction is used in the form suggested in Ref. [pq| . The point-like potential 
in the configuration space would lead to the divergences. For this reason, the following 
regularization of the potential is considered 



Vr(r) -> AgV T ,reg(r; A) , V T>reg {r; A) 



(AV5F) 
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'A 2 



In the momentum space, we have 



Vr,reg(p - p' , A) = exp 



A 2 (p-p') 2 

4: 



Now, using the following representation of the Legendre polynomials 

1 d n 



Pjz) 



2 n n\ dz 1 
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with the use of the identity 



e az z n dz 
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after the partial-wave expansion of the t'Hooft potential we obtain 

J^^reafop'; A ) = S LoVT, reg (p,p'] 0) , 

that reflects the point-like character of the t'Hooft interaction. Here, 

A 2 (P 2 +P /2 ) 
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(4.31) 



(4.32) 



(4.33) 
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(4.36) 
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In accordance with the Eq. (|4.35|) , all partial waves except L = in the partial- wave expan- 
sion of the t'Hooft potential are neglected even at nonzero A. 

As it was mentioned above, the t'Hooft interaction was introduced in order to provide 
the mass splitting between the pseudoscalar and vector octets within the framework of 
the constituent quark model, which in QCD is due to the spontaneous breaking of chiral 
symmetry. The quantity g that appears in Eq. (|4.31| ), is the matrix in the flavor space. The 
matrix elements of this matrix between various meson states 



it, p, K, lo = rjr, 



V2 



[uu + dd) , $ = rj s = ss . 



(4.37) 
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(ir\g\ir) = -g = (p\g\p) , (K\g\K) = —g' , 

(Vs\g\Vs) = , (Vn\g\Vn) = g , (Vn\g\Vs) = V2g'{rj s \g\r] n ) 



(4.38) 



where g and g' are two independent coupling constants that are considered to be the free 
parameters of the model. 

The 7] and rj mesons are the superpositions of r\ n and r] s . In order to take the mixing 
into account, we introduce the matrix notations 



$M S (P) 



Va(p,p') 



Mp,p') 



G (M B ;p) 







fn,M s (p) 
$ S ,M B (P) 

K,a(p, p') 

K,a(p,P' 

(p,p') Ks,t(p, p') 

K«,t(p, p') V aa ,T(p, v ) 



(G nfi (M B]V ) 

I G s , (M B ;p) 



A = C,OG, 



(4.39) 



The radial wave functions R^qqq (p), f = n,s describing the 77 and rf mesons, obey the 
following system of equations 



[M B - (a lWl + a 2 w 2 )} f Rf^\p) = A { p a * ] {M B] p) ]T / p' 2 dp' 



x 



, , J 



(A/^ 2) (p)^ 
/' L 

+ ai a 2 U^ a2 \p) + ^a' 2 Af^t a ' 2 \p'WteM^ V(Vf\mVf')W£$V)} ■ (4-40) 



The functions R^*^\p), f = n,s satisfy the normalization condition 
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p 2 dp 
(2^)3 



E 



&2 2) (m Bj p)(r^\p)) + &^\m b ,p)(r^\p) 



2M B , (4.41) 



n, s 



where M B is either M n or M v >. 

The equations for other mesonic states can be obtained, replacing (r]f\Ag\T]f), f,f 
by the corresponding matrix elements from Eq. (|4.38|) . 

Note that the mixing in $ — uj and 77 — 77' systems has been recently also investigated in 
Refs. [[41]] within the Nambu-Jona-Lasinio (NJL) model, with an account of the relativistic 
confinement potential (Lorentz vector structure only) and the t'Hooft interaction. 



D. Solution of the equations 

One has to specify the numerical procedure for the solution of the system of radial 
equations ( |4.16| )-( fLl7| ). A possible algorithm looks as follows. One chooses the known 
basis functions denoted by R n L(p)- The radial wave functions are expanded in the linear 
combinations of the basis functions 



3 rf(aia 2 )/ 



n=0 



„(aia 2 ) 
-nLSJ 



RnL(p) , 



(4.42) 



where ' 2 } are the coefficients of the expansion. The integral equation for the radial wave 
functions is then transformed into the system of linear equations for these coefficients. If the 
truncation is carried out, the finite system of equations is obtained that can be solved by 
using conventional numerical methods. The convergence of the whole procedure, with more 
terms taken into account in the expansion ( 4.42|) , depends on the successful choice of the 
basis. In Refs. P4]P3| , |25| - P^ , P^ ] , where the linear confining potential is assumed, the basis 
functions are chosen in the following manner 



Rudy) = N nL y L L 2 n L+2 (y)e 



-1//2 



y = PP: 



(4.43) 



where L 2 ^ +2 {y) are the Laguerre polynomials, and /3 is the free parameter. In Refs. []I3 , |27 , 55 
the non- relativistic oscillator wave functions (again containing the free parameter), were used 



in spite of the fact that the linear confining potential was assumed. In Refs. [p4],p!7Hl9|,pl|,|35 



the same basis functions were used, but without the free parameter, due to the fact that 
the confining potential was taken in the harmonic form, with the parameters already fixed. 
Finally, in Ref. ||26|| , the harmonic oscillator basis was used, whereas the confining potential 
had the general form given by Eq. ([OO). 

To clarify the choice of the basis functions, let us consider the non-relativistic limit of 
the equations (|3.61|) . In this limit, one can replace 70 — > 1, 75 — > 0, ■y — > 0. Consequently, 



V^Vog + V c + V T 



(4.44) 



Further, to derive the non-relativistic limit of the equations, we expand the kinetic term 
aiWi + a2W2 in Eq. ( |3.61| ), retaining terms up to and including 0(p 2 /m 2 ). In the right-hand 
side of this equation, the function A^ aia2S> {M B ]p) can be replaced by its value at p — 0. In 
a result, we obtain 
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: A( ++ )(M B ;0) = 1, 



(4.46) 



CJ : A<++)(Mb;0) = 1(i + ^). 

2 V m 12 / 



(4.47) 



The non-relativistic limit in the MNK version is more tricky. For a general Mb, there 
emerges an arbitrary function of the ratio Mb/tti^. However, if one uses the non-relativistic 
approximation also for the bound-state mass Mb = w,\2, then A^ ++ \Mb] 0) = 1 for both 
the CJ and MNK versions. Below, we shall use this approximation. 
Since in the non-relativistic limit (see Eq. ( |4.3| )) 

*£#*(P) = X£#*(P) , (4-48) 



the non-relativistic limit of Eq. (|3.61|) with the harmonic confinement potential (|4.22|) only, 
is given by 



(4.49) 



Performing the partial- wave expansion of Eq. ( [4.49 ), we obtain the equation for the radial 
wave functions 
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Rl(z) = 0, (4.50) 



where 2; = p/p, and p = y /in^oy |«s( m i2)- The solutions of this equation with the energy 
spectrum 



e B ] = ~ots(ml 2 )V + J^a s (m 2 12 )uj (2n + L + ^) , 



(4.51) 



are the well-known harmonic oscillator wave functions 

( z 2 \ 3 
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(4.52) 
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where \F\ denotes the confluent hypergeometric function. 

The functions R n L(p) can be used as a basis for the expansion in the general case Q4.42|) . 
The system of equations for the coefficients is given by 

M B & = E EE , (4-53) 

where the matrix H(Mb) is given by the convolution of the potential and various kine- 
matic factors that appear in Eqs. ( |4.16| )-( |4.17| ), with the wave functions of the basis. From 



Eq. ( 4.53 ) it is immediately seen that, in general, the eigenvalue equation for Mb is not a 
linear one, and should be solved, e.g., by iterations. 

In order to actually solve the system of equations (|4.53|) , one has to truncate it at some 
fixed n = N max . then, c^jn are determined from the system of A(N max + 1) (2(N max + 1) 
in Salpeter and Gross versions) linear equations. This procedure determines the eigenvalue 
Mb as well, either directly, when the matrix H(Mb) does not depend on Mb, or by using the 
iterative procedure. Having solved the eigenvalue problem at a fixed value of N max , one has 
then to check the stability with respect to the change of N max - if the calculated eigenvalues 
do not converge with the increase of N max , the original system of integral equations is 
declared to have no solutions. 

Note that the system of equations (|4.53| ) is homogeneous in c^gj^ . This means that the 



solution of the eigenvalue problem determines these coefficients up to an overall factor that 
can be fixed from the normalization condition 

E Efe } fe } rp 2 d P fir 2 \M B ;p)R n L(p)Rn>L(p) = 1. (4.54) 
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In the CJ and MNK versions, the function f[ 2 \Mb',p) has the second-order pole, so in the 
normalization condition one encounters singular integrals of the following type 



f(x)dx 
o (x — x ) 2 



'(so) = L , (4-55) 



where f(x) is the regular function that obeys the conditions /(0) = f(oo) = 0. The integral 
in ( 4.55|) can be regularized according to 



f{x)dx _ /-2-o (f( x ) - f'{x ))dx + r°° f{x)dx ^ 



(X — Xq) 2 JO X — X J2x x — x 

The first question which one may be willing to investigate, is the manifestation of the 
Lorentz structure of the confining interaction in the bound-state mass spectrum, especially 
in the case of light quarks. This question was addressed, e.g. in Ref. [Q, where the scalar, 
timelike vector, and their equal- weight mixture were studied on the basis of Salpeter equation 
(this corresponds to the choice x = 0; 1; 0.5 in Eq. (|4.2j ), respectively). It was demonstrated 
that the stable solutions to the Salpeter equation in the light quark sector do not exist 
for the scalar confining potential x — 0, and do exist for x = 0.5 and x = 1. Further, in 
Ref. [0 , the structure 7W g) yH 2 ) was considered as well - it was demonstrated that in the 



case the stable solutions do not exist. In Ref. [19], more general conclusion was obtained 



it was demonstrated that the stable solutions in the light quark sector exist for any x from 
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the interval 0.5 < x < 1. This result was confirmed in Refs. p7|j31[| . Further, in Ref. |19|], it 
was shown that in the heavy quark sector nothing really depends on the mixing parameter 
x - the solutions exist everywhere and practically do not change when x varies in the whole 
interval < x < 1. This result is easy to understand. Indeed, the projection operator 
(A^^ — A^ 2 jTo^ ® 7cj 2 \ that is present in the Salpeter equation, in the heavy quark limit 
is equal to |(7o + 1^)1^ ® so that the confining interaction in this limit equals to 

(AS+ +) " Afe- ) ) 7 f ) ® 7S 2) [x7S 1} ® 7o W + (1 " ® I (2) ]V C (P ~ P) 



-^(7S 1) +7S 2) )^c(p-P / ) (4.57) 

at mi, r«2 — ► oo, and does not depend on a; at all. Note that in the literature we encounter 
the different choice of the parameter x: x = 1 x = °- 5 HSH, z = ©H,ll 



Note also, that, as it was shown in Ref. f2l| , the non-existence of the stable solutions at small 
x in the light-quark sector is related to the presence of the "negative-energy" component in 
the Salpeter wave function. 

The same question can be studied in other - GR, MW, CJ and MNK - versions, that, 
unlike the Salpeter equation, have the correct one-body limit. For the MW and CJ versions 



the investigations were carried out in Ref. f23fl . Here, the problem was studied in the 



configuration space, and for the confining potential the following Lorentz structure was 
assumed: Vc(p,p') = [xj^ ® Y {2) + (1 - ® I {2) ] Vc(p - p') where for V" c (p - p') 

a linear form was chosen. It was demonstrated that this potential should be "more scalar 
than vector" in order to provide the existence of the stable solutions. More detailed study of 
MW, CJ and MNK versions in the momentum space was carried out in Refs. |Tl|j35|] , where 
the harmonic confining potential was used, with the Lorentz structure given by Eq. (|4.2| ). 
The following states ds : %, 3 Si, 1 Pi, 3 Po, 3 Pi, 3 P2, 1 P2, 3 P>i, 3 D 3 , cu and cs : ^i, 1 Pi, 3 P 2 , 
were considered. It was demonstrated, that in all versions the solutions always exist at 
x = 0, whereas for x — 1 for the majority of the states there is no solution. This is just the 
opposite to the Salpeter equation case (see above) - there, at x — 1 there are the solutions, 
whereas at x — 0, the solutions for majority of states cease to exist. Put differently, the 
existence/nonexistence of the solutions depends critically on the value of x, and the criteria 
vary from version to version. In addition, the criteria depend on the details of the potential 
- in particular, the strength parameter u introduced in Eq. (|4.22|) . Note that the instability 
mentioned, is now caused by the admixture of the mixed (H — ), ( — h) frequency components 
in the bound-state wave function. One may look for the admissible window in the parameter 
space, where the solutions of all versions simultaneously exist and approximately coincide. In 
this way, one may judge on the Lorentz structure - assuming that the whole physical picture 
of the qq bound states based on the 3D reduction of the BS equation, is viable. From this 
study, one has to reject the MW version that poorly agrees either with other versions or 
with data. Further, on the basis of SAL, CJ and MNK versions, one can determine the 
acceptable interval for the mixing parameter x: 0.3 < x < 0.6. 

Both the fine structure (P-wave splitting), and the hyperfine structure ( 3 S'i — 3 D\ split- 
ting) of the qq states depends on the value of the mixing parameter x. As it was shown in 
Refs. p5| , p3| on the basis of Salpeter equation, the spin-orbit splitting in the light quarkonia 
can only be described by the mixture of scalar and timelike vector confinement. However, 
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as was shown in Ref. ||33|| , the fine structure and the hyperfine structure can not be si- 
multaneously described by simply varying the value of the mixing parameter. Finally, in 
Ref. ||35|| , more general - and pessimistic - conclusion was drawn: neither of the versions - 
SAL, MW, CJ or MNK - with the dynamical input specified above, does not describe even 
qualitative features of the whole mass spectrum of qq bound states with x inside the interval 
0.3 < x < 0.6. Clearly, the problem calls for the further investigation. 



V. DECAYS OF THE MESONS IN THE CM FRAME 

Further information about the bound qq systems may be gained, investigating their 
decays. Below, we consider exclusively the decays that proceed into the CM frame of the 
bound state?]. These are: the weak decays of the pseudoscalar mesons P — > \iv, the leptonic 
decays of the neutral vector mesons V — > + e~, and the two-photon decays M — > 77. The 
corresponding characteristics are: the weak decay constant fp, the leptonic decay width 
T(V — > e + e _ ) (or, the leptonic constant fv), and the two-photon decay width T(M — > 77). 

The expressions for the quantities fp and T(V — > e + e~) were obtained in Refs. 
T8| , |25|p8| , |29 | in the framework of Salpeter equation, directly in terms of ^(^(p) = $ aa (p) ± 



$6&(p), or ^ aa (p) = $ (++) (p), ^ bb (p) = <& ( ~ _) (p) (see above). In Ref. these quanti- 
ties were evaluated in the framework of SAL, CJ and MNK versions written in the form 
( |4.5|) -( [i~l~8l) , that corresponds to the representation of the wave function in the form ( |4.3|) - 



( ff.4| ). The main conclusion that comes from this investigation, is that the results do not 
depend much on the choice of the different 3D reduction scheme. The quantity T(M — * 77) 
was evaluated in Refs. PB|J2D|J3^J3^1 for the systems (71,77,77'). Below, we shall follow the 
derivation presented in Ref. ||35|| . 

For the calculation of the quantities listed above, we need the wave functions ^lsjMj (p) 
which, according to Eq. (fO|), is expressed via Xl5jm 7 (p)- The partial- wave expansion for 
the components of the wave function reads 

= (n\LSJMj)Mr 2 \p) R ( ls?\p) , 

[®LSJMj(P)\ab = ~( Sn - CTn ) HLSJMj) R> LS j '{p) , 

W2 -T OC2 m 2 

./■(0102)/ \ 

r^(«i a 2) / M / \ / rcTi / \ a lP J *l2 \P) j D (a 1 a 2 ) , \ 

[®LSJMj(P)\ba = ( sn + CTn ) (n\LSJMj) Rlsj '{p) , 

r*(aioa) / m s i2 + v {1) (T {2) , lTQT , /T \ a l a2P 2 N\ L 2 ia2 \p) ( ttia2 ), > , . 

= 3 (n\LSJ Mj ) {wi + aimi){w2 + a2m2) Risj \P) , (5-1) 



2 The treatment of the decays which can not be confined to the CM frame, implies the specification 
of the Lorentz-transformation rules for the instantaneous potentials and 3D wave functions. Due to 
the Lorentz covariance, the dependence on the 0-th component of the relative momentum emerges 
into the transformed wave functions, that renders the problem extremely complicated, and the 
further assumptions are necessary. We do not consider such processes here. 
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where S = |(<r« + <x( 2 )), a = |(<r« 
by Eq. ($T§. 

Using now the identity which is valid for any operator O 



cr^ 2 )), and the operators S12 and er^cr^ are given 



6(0- 



(1) 



n)(n\LSJMj) = £ (n\L'S'J'Mjr) (L 1 S 1 J 1 Mj>\Q\LS JMj) , (5.2) 



L'S'J'Mj, 



and the well-known expressions for the matrix elements of the operators Sn, <xn, S12, from 
Eq. ( |5.1| ) it is straightforward to obtain 

(n\LSJMj)Mt a2) (p)R ( LSJ 2) (p), 



LSJAfjVPJJaa 



T 1 



J+ 1 
J 



J 
J + l 



N 2J+1 (n| J + UJMj> + N 27Ti (n| J - 11JMj> 



X 



[*S±Ti 2 k(p)] a6 



j 

j + i 



j + i 
j 



N 27Ti (n| JUMj) T N 2JT1 (n| J0JMj) 



X 



J + l 
J 



J 

J+l 



■\ 2J+1 ( "l J + 11 JA « ± \ 17TT <n| J - 11 JM ^ 



X 



T,(«l"2) / NT 

J+llJMjVPJJ&a 



J 
J+l 



J+l 
J 



66 



N 2J+1 J1JMj ) T \ 2J+1 ^ n l J0JMj ) 



x 



J±llJMjlPJJ66 



"2J+ 1 



2JJ(J+1) 

(n| J ± HJMj) - ) (n| J + llJMj) 



2J+ 1 



x 



xaia > A4- ai - (, ' ) (p)jfej(p). 



With the use of these expressions, we can explicitly calculate the quantity 



. 2 ("102) 



V) 1 



i,j = a,b, 



(5.3) 



(5.4) 



as an 2 x 2 matrix in the fermion spin space. To this end, we explicitly introduce the fermion 
spin coordinates (j\ and cr 2 (cr, = ±|, z = 1, 2). Then, we have 
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(n\LSJMj) = {no l( j 2 \LSJMj) = ]T (LSm L m s \JMj) (^a^Sms (n\Lm L ) . (5.5) 

m L m s 

Further, with the account of the following relations 

1 



<^ 2 |00>H^) = i(J \ 



I = 00, 



(^-(7 1 (7 2 |l0)(-i(7 2/ ) 
.11 

^2 2 " 10 " 2 ' 1 ± ^("^ 



J_ /l 

V2V0 -u - ^ 



(V 

(?) o 



1 ( ■ ^ 1 

2 (To - * - wt„; = CT± = </?l±l , 



it follows that 

(nlLSJMj) = J2 (LSM L m s \JMj) (n\Lm L ) Sms = «n|L) <g> ^ s ) JMj • 



m L m s 



(«i«2) 



For the quantity $l 5 jMj(p) = E Ql a 2 *l5JMj(p) we obtain 



(5.6) 



(5.7) 



[*,(») jm»L = [=F)| ^7_f «n| J + 1) ® ft)™' + «n| J - 1) ® fc) 



J + l 
J 



X 



J±llJMjVP)] a & 



X 



E«2M ( 2 1_Q2) (P)^ 2 ](P), 



j 

x j + i 

\ 2J+1 



((n|J) (g)^ 



/ J+: 
\ 2J+1 



((n|J)®^ ) JMj 



x 



[*j(°)jM>)] 



6a 



/J+l 



J 
J+l 



X 



\| «n| J + 1) ® ^ ± ((n| J - 1) ® fc) 

E«4 aia2) (ri^(i 



X 



[^J±llJAf J (p)] 6o = 



J 

J+l 



J+l 
J 



X 



\ 2lTl ((n| J) 0l)JMj T \ 2J+1 ((n| J) ® 0o) 



JMj 



X 
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[$J±UJMj(p)} bb 



± (n|J) 



<Pi 



E ^4 ai " a2) (p)tt(p) 

QlCt2 ^ 1 ' 



— «n| J ± 1) ® ^) JMj - 2 V J(J + 1 ) (( n | j zp i) 



x 



2J + 1 vv 1 ' T ' 2J + 1 

x £ a^X^^^P)^! 2 ]^) • (5-8) 



In order to evaluate the constants fp and fy, we need the bound-state wave function of the 
qq state at r = 

^LSJMj{r = 0) = ^L5JA/. 7 (r = 0, C7i, 0- 2 ) = J *LSJAfj(p) = j J^y ®LSJMj(p, <T\<T 2 ) , 

(5.9) 

where, according to Eq. (|2.26|) , 

vf/ ( n \ — { (®LSJMj{p))ab {^LSJMj{p))aa\ /r im 

^LSJMj{P) ~ / f / n (x y ( I • 15.1UJ 

The decay constants /p and /y for the pseudoscalar (L = S = J = 0) and vector 
(L = 0, S 1 = J = 1) mesons, respectively, are given by 



S^Mgfp = \/3 tr [^ 0000 (r = 0) 7 "(1 - 7 5 )] , 

fvW = v / 3tr[§oiiA(r = 0) 7 At K(A), A = ±1,0, (5.11) 

where the factor a/3 stems from the color part of the wave function, and e M (A) is the 
polarization vector of the vector meson P7| 

e M (A = ±l) = T-7=(0,l,±i,0), e M (A = 0) = (0,0,0,1) , in CM frame. (5.12) 
V2 

Now, using the equations ( |5.9|) -( [5TTT1) , we obtain 

fp = ^ r E [A/g^GO - W/^ 2 ^)] 4oo Q2) W , (5-13) 

Mp Jo (2vr) d aiQ2 

/v(A) = -c^ov^ r E [^V) + ^^""^(p)] 4ir 2) b) = <Wv • 

Jo (27r) J QiQ2 3 

The leptonic decay width of the vector mesons (p°,uj, $) is given by 

-, (5.14) 



where 



Ml 3 A= tl i0 '^ v 71 3M| 



.2-2 - „ / „ -2 ,111 



e 9 /e, e 9 = (-, — , -) , (5.15) 
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for p°, to, $ mesons, respectively. Here, e q denotes the expectation value of the quark charge 
in the units of the elementary charge e. 

In order to explain the reason, why the quantity e q appears in the expression ( |5.14 ), let 
us note that the leptonic decay of the vector meson in the lowest order in e is described 
by the diagram depicted in Fig. 1. Taking into account the flavor structure of the wave 
functions 



p° ~ — y= (uu — dd) 



uj ~ —j= (uu + dd) 
V2 



$ ~ ss . 



(5.16) 



we obtain, that the transition amplitudes of the vector mesons into the photon are propor- 
tional to 



(u -> 7) 



3V2' 



($-7) 



(5.17) 



from which the Eq. ( 5.15Q follows directly. 

Further, taking into account Eqs. (|4.42|) , ( |4.52|) and (|4.54| ), we can express the quantity 
fp and the leptonic decay width in terms of the dimensionless wave functions R^gj 2 (z) 



fi 



T(V 



V6f /2 



B 



E 



z dz 



E 



R 



(aia2) , 
000 I 



(5.18) 



Ot\OL2 



z dz 



Mfr a2 \p,z)+^Mr i - a2) (P,z) 



RiT 2 \z) 



n=0 



where the functions 

RtsT\z) 

satisfy the normalization condition 

E / z 2 dzffr a2) (M B ;p,z) 



E ^LSJn R nL(z) 



(5.19) 



R { LSJ 2 \z) 



1 . 



(5.20) 



Next, we consider the two-photon decays of the neutral mesons. The amplitude of the 
two-photon decay of the qq bound state with equal-mass quarks in the lowest order in the 
coupling constant e is given by the diagrams depicted in Fig. 2. Th the CM frame where 



P = (M B , 0), this amplitude is equal to p4 

d 4 p 



T(X 1 X 2 



(27T) 



(5.21) 



where ^ m b (p) is the BS amplitude of the qq bound state which satisfies Eq. ( |2.25| ) and is 
written in the form (|2.26j) . Further, k\ = k) and k 2 = (^f 2 -, — k), where k is the relative 
3-momentum of photons in the CM frame directed along the z-axis, and 6j = e(A») are the 
polarization vectors for the emitted photons. Due to the fact that the emitted physical 
photons are transversely polarized, one needs to consider only the values \ = ±1 for which 
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e(A! = ±1) = tA= (0, 1, ±i, 0) , e(A 2 = ±1) = ±^= (0, 1, T i, 0) , (5.22) 

and 

AM = ±x ) = 7f + ^) ' ^ = ±1} = 7f (T7 * + i7 «) • (5 - 23) 

Further, one may rewrite the expression, entering the integrand in Eq. (|5.21|) in the following 
manner (Below, we follow the derivation given in Ref. |f48|.) 

f M (p, k) = -i ^(f +p-h)6-i &S(^ + P -k 2 ) i x (5.24) 
a«(p-k) | ai 2 ) (p-k) | 4l } (p + k) | 4 } (p + k) 



p -w(p-k)+i0 p + w(p - k) - iO p — w(p + k) + iO p + w(p + k) - iO ' 



where w(p ± k) = ym 2 + (p ± k) 2 , and 

a$(p - k) =^A (Q) (P - k) 7 o f6 , 4i } (p + k) = / 2 A (a) (P + k) 7o h ■ (5.25) 

Note that, of course, the relation of the BS amplitude ^ m b (p) an d the 3D amplitude ^m b (p) 
is different in different versions of the 3D reduction. In particular, in the Salpeter version, 

= S(~ +p)T{p)S(~ +p) , (5.26) 
where, taking into account Eq. fl2.26|) , we have 

On the other hand, from Eq. (|5.26|) one may obtain 

r(+-)( P ) / i i 



®m b (p) 



+ 



M B - 2w + iO \p - ^f- + w - iO p + ^ - w + iO 

r(-+)( P ) / 1 1 



M B + 2w \p - Ma + w _ iQ po + Ma _ w + jQ 



+ r<++)( P ) ^ i i 



where 



Mb Vpo - - w + i0 Po + ^-w + iO 

+ ^M(^1 ^ ), (5.28) 

Mb \p -^ + w-i0 p + ^ + w + iOJ 



r M) (p) = A W(p) 7o r(p) 7o A^(-p) . (5.29) 
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After integrating Eq. ( |5.28|) over p , we obtain the Salpeter equation for the equal-time 
amplitude 



*M B (P) 



zr(+-)(p) ir(-+)(p) 



M B -2w + iO M B + 2w 



(5.30) 



Now, substituting ( |5.28|) into the expression of the two-photon decay amplitude (|5.21|) and 
integrating over p , we obtain: T(±=)=) = and 



T(±±) = iV3e 2 q J 



1 



d 3 p 

(2^)3 \ \M% -(w + w(p-k))i 



tr 



+ fl( r (++)(p)_ r (--)(p))+ , 



r ( - +) ( P ) + n 



i 

2 



(r( ++ )(p)-r(-)( P ))( 7oT757 ,) 

+ w + w(p-k))W /s (p)) x 



X 



u>(p — k) 



((1 ± 75707*)*" + (jz =F 757o)(p 2 - fc)) 



iMl-( W + W (p + k)) 2 



tr 



(r( ++ )(p)-r(-)( P ))( 7 oT757, 



+ U( r(++) ^ +r( " ) ( p )) + VM B + 2. 



Mf 



r ( - +) ( P ) + K^f- + w + w( P + k)))^ Ms ( P )) x 



i 



X 



w(p + k) 



((1 =F 757o7^)^ + (7* =F 757o) ip z + k)) 



(5.31) 



For the further transformation of this expression, one may use the fact that the static 
potential V(p, p') has the Lorentz structure given by Eq. ( |4.2j ). Then, 



V(p, p , )^M s (p') = V OG (p - p')7o^m s (p' )7o + Vc(p - p') f xry$ MB (p')7o 
+ (1 - x)V MB (p')) + V T (p - p')4£((tr *M B (p')) + 7str (tt^p'frs)) . 



(5.32) 



From this, one can directly obtain 

d 3 p' 



r (^)( p ) = -za( q )( P ) y 
+ vb(p-p'; 

+ v T (p - p')4g 

where 



(2*0 



Vog(p - P 



/^ I &ab(p') 0aa(p') 
W) L(P') 



ab (p') (2x - l)0 aa (p') 

(2x-l)^(p') i ba (p') 

tr (<L(p') + <L(p')) -tr (<L(p') + 66 (p')) 
"tr (L(P') + <L(p')) tr (0 afe (p') + 4> ba (p>)) 



(5.33) 
AW(-p), 



1/3 (PJ = -«0a/KPK 



(5.34) 
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are the components of the meson amplitude in the spin space. After performing the partial- 
wave decomposition of these amplitudes, the expression for the quantity T(\ 1 \2) takes the 
form (note that we have replaced t'Hooft interaction by its regularized version). 



1 S'o state: 



T(±±) = ±^L [ 
y ! (V2ti) 5 Jo 



(V27TY 



p 2 dp 



■ur 



p dp 



2mw T . , , . 2mw r . , , 
M B (M B + 2 a ) J °( M ^ + M^ P h(MB ^ 



2mw 



I 2 (M B ;p) )V^p,p')R 000(ab+ba) (p') - (2x -!)((-£- + 



m 



M B p(M B + 2w) 



)J (M B ;p) 



2w 2w 2 \ 

+—I 1 (M B ;p) + j^I 2 (M B ;p)JvS{p,p')B<m(aa-bb){p') 



+(2x-l) 



mw 



+4g 
Aw 2 



p(M B + 2w) 
2mw 



J (M B ;p) Vg(p,p')Rooo(aa+bb)(p') 



+ 



x = l,Vc^V OG 



p(M B + 2w) 



J (M B ;p)V^ reg (p,p';A) R 00 o(aa- bb) (p' ) 



+ 



m T , , „ , s . Mb t0 



M b 
3 P state: 



— h(M B ;p)Rooo(ab+ba) (p) + (^— h ( m b~, P) + h(M B ;p))R 00 o(aa-bb) (p) 
p ip 



T(±±) = ± ^L rt** 

y 1 (V2n) 5 Jo w 2 



p dp 



2m , m 2 



+ 



P 



M B pM B M B + 2w 



)Jo{M B -p) 



MBp' ir "" r ' ' M B (M B + 2w) • M 2 B 



(5.35) 



+^/ 1 °(M B ; P ) + _ T |^L_/ 2 (M fl ;p) + ^I 2 2 (M B ;p) )V»{p,p')R mab+ba) {p') 



•(2s-l)( 



4m 2 w 



-J (M B ;p) + —I°(M B ;p) + (— + — — -)I 2 (M B ;p) 

M 2 (M B + 2w) V ,FJ M B 1V Mb p(M B + 2u>) 



m 



+2(£ 2 



M| M B (M B + 2w) 



/|(M s ;p) )^ c 1 (p,p / ) J R 110(aa _ 56) (p / ) - (2x - 1) x 



(J (M B ;p) - J 2 2 (M B ;p)) - h{M B -p) ) ^(p^O^noCaa+wjO/) 



M s + 2w M| 



+ 



x = l,V c ^V c 



OG 



+ 







m 



m bP ^ m ^ + mIv MMb '' p) ' 2M B 

"'" Jo(M B ; p)-H- h(M B ;p) - /|(M B ;j>) 



m 



h{M B -p) 



+^l!(M B ;p)) + 



M B + 2w\M 2 " u ^'^^ J Mp' Ai " fl,J,y M 2 



x 



4w 2 



X^T.res^P; A )^110(a6+6a) ) ) +T7 ~ *3 ( M B! p) RllO(ab+ba) (p) 



M 



B L p 



+(^-I 3 (M B ;p) - l!(M B -p))R U0(aa _ bb) (p) 



(5.36) 
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where 

RLSJ(aa+bb)(p) = R { LSJ(P) + RjLSj' ' (?) , 

RLSJ(aa- bb) (p) = -{R% + J\P) - R { LSKP)) + "(ifefo) + 4^(p)) , 
i?X5J(a 6+6 a)(p) = ^(4^(P) - + ^(4 + ^(p) + R { LS + ;\P)) , (5-37) 



and 



j . 2{w + w + )-M B 2{w + w-) - M E 
Jo = m — : — In 



J = ln 



2(w + «;+) + M B 2(w + wJ) + M B 

2w(w + w + ) + M B p 
2w(w + u>_) — M B p ' 



T - 1 T W T 

Aw w 2 _ w T 

u; + + W- M B p 2p 

2M B w 



w + + w_ p 

2 _ 2w 2 5w 2 -\Ml + w + w- w 2 
h — I" q 7 : \ I T Jo , 

p 3 P{W + + U>_) p z 

1° = ln 2 ( W + W +)~ M B 



2{w + w^) -M B 
v 

W + + W- p 



M B - 2w w 
1 H : — ^3 j 



, 2 _ w (M B - 2w)(w 2 + \M 2 B - w + w. + 3M B w) w 2 
p SMbP{w + + W-) p z 



with w± = Jw 2 + \M 2 B ± MbP- 

Now, we consider other versions of the 3D equations. Since we consider the equal-mass 
case, the Gross equation can not be used. For this reason, we shall restrict ourselves to 
study of two-photon decay processes in CJ and MNK versions. In these versions, there 
exists a relation between 4D and 3D free Green functions given by Eqs. ( 3.46Q and ( p. 54 ). 



This relation can be immediately translated into the relation between the 4D and 3D wave 
functions 

#i^(p) = 27ri$(po)*M B (p), (5.39) 

where for the MNK version the equality p$ = holds for the equal-mass case. As to the MW, 
version, here the relation between \l/ Mg (p) and $m b (p) does not exist due to the definition 
of G^ w (Mb,p) (|3.31|) . For the above reasons, below we restrict ourselves to the CJ and 
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MNK versions only. Substituting the expression ( |5.39|) into Eq. ( |5.21|) , with an account of 
G5.24j ) one obtains 



r(AiA 2 ) = iVZe 2 q J -^tr|f Mfl (p) 



(2tt)' 



4+ } (p-k) a^(p-k) 
tu(p — k) w(p — k) 

a£>(p + k) afe^p + k) 



(-), 



w(p + k) w(p + ~k) 



From this equation one readily obtains 
T(±,±)=±iV3e 2 q J ^tr{* MB (p 



(1 ± 7 5 7o7z) m + {lz =F 757o) (g* ~ fc ) 
u> 2 (p — k) 

(1 =F 757o7^)™ + {lz ± 757o)(p* + k) 



w 2 (p + k) 



(5.40) 



(5.41) 



Substituting ^ MB (p) in the matrix form given by Eq. ( |5.10|) , we finally obtain for the CJ 
and MNK versions 



1 S state: 



2 r°° 



pdp 



-2 



w 2 + \M% 



M B p J ( M B^))j^Rooo (ab+ba) 



rn 



(P) 



+ 



fly w 2 - \M 2 B 



\M. 



B 



2M 2 



J{M B -p)\R 



B 



£ 000( aa—bb) 



(P) 



(5.42) 



5 P state: 



T(±±) = -zv^e 2 , 



(2vr) 5 / 2 Jo 



+2 



^ ( M B I P) T7~ R H0(ab+ba) (?) 

M B 



w 2 - \M 2 B w A - hM% - 



M 2 



2 M | p J( M B;P) )RllO(aa-bb) (P) 



where 



J{M B ;p) =ln 



w 2 + \M% + M B p 
w 2 + \M\ - M B p ' 



(5.43) 



(5.44) 



It is important to note that in the Salpeter version the two-photon decay amplitude depends 
on the potential both directly and indirectly, through the radial wave functions, whereas in 
CJ and MNK versions this dependence enters only through the radial wave functions. 
For a given meson, the two-photon decay amplitude can be rewritten as 



T(AxA 2 ) = e 2 eJ je// v / 3T(A 1 A 2 ; LSJMj) . 
The decay width is given by 



(5.45) 
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a 2 1 



r(meson -> 77) = 3tt — 51 g J ;/ f (AiA 2 ; LSJM, 

M B 2(2 J + 1) AiAaMj 



./J , (5.46) 

where e q e ^ depends on the choice of the meson flavor wave function. If this function has a 
simple form qq then e q e ^ = e q . However, if the meson wave function is made up of different 
flavor states aq\q\ + /3q2<l2, the expression for *? qe t* is more complicated. Consider as an 
example calculation of this factor for 7r° and r\ n states. The flavor structure of the wave 
functions is given by 

7T° ~ —j= (uu — dd) , r] n ~ —j= {uu + dd) . (5-47) 
V 2 V 2 

It follows then straightforwardly that e qe jj- = ^= and e q>e ff = for tt° and r\ n states, 
respectively. Further, the decay amplitudes for the physical r\ and 7/ mesons are the linear 
superposition of the ones corresponding to r\ n and r\ s ~ ss states. 

Note that the two-photon decays of 7r°, rj, rj' mesons was also studied in the NJL model, 



taking into account the relativistic confinement and the t'Hooft interaction [42 . 
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FIGURES 




FIG. 1. Decay of the meson into electron-positron pair 




